Crossover from one to three dimensions for a gas of hard-core bosons 
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We develop a variational theory of the crossover from the one-dimensional (ID) regime to the 3D 
regime for ultra-cold Bose gases in thin waveguides. Within the ID regime we map out the parameter 
space for fermionization, which may span the full ID regime for suitable transverse confinement. 
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It has recently been proved by Lieb and Seiringer |l| 
that in a suitably-defined dilute limit, the many-body 
ground state of a trapped ultra-cold gas of bosons in two 
or three dimensions exhibits Bose-Einstein condensation 
into that orbital which minimizes the Gross-Pitaevskii 
(GP) energy functional, and in fact that the condensa- 
tion is complete in the sense that the condensed frac- 
tion is unity. Their work should be consulted for precise 
definitions and hypotheses required for the proof; here 
we reiterate only a few points relevant here. In the 3D 
case the dilute limit is defined as a — > and N — > oo 
with both the trap potential and Na fixed where a is 
the s-wave scattering length and N the number of par- 
ticles, and in the 2D case it is defined as a — > and 
N — > oo with the trap potential and N/ \ \n(a 2 N)\ fixed. 
They point out that in ID their proof fails, where there 
is presumably no Bose-Einstein condensate (BEC) even 
at zero temperature (many-body ground state). In fact, 
the exact many-body ground state of a spatially uniform 
(untrapped) ID Bose gas with repulsive zero-range (delta 
function) interaction was found long ago by Lieb and Lin- 
iger (LL) [Q] for all values of a ID coupling constant giu, 
and was shown by them to reduce for giu — > oo to the 
exact many-body ground state of the impenetrable point 
Bose gas ( "Tonks gas" ) found previously by one of us jjj , 
for which Lenard proved rigorously |Q that the occupa- 
tion of the lowest orbital is bounded above by const.-^N, 
ruling out BEC (occupation proportional to N). In the 
case of a trapped Tonks gas no such rigorous bound is 
known, but our numerical evaluation of the largest eigen- 
value of the reduced single-particle density matrix of our 
exact many-body ground state Q suggests strongly that 
the most highly occupied orbital has occupation behav- 
ing like N p with < p < 1, again indicating absence of 
true BEC. In what follows we will use the term ID con- 
densate to describe the ground state of the trapped gas 
when the the system is still ID, in that it has the trans- 
verse profile of the trap ground state, but the energy has 
deviated below that for an impenetrable Tonks gas. 

It is clear from the above discussion that for real atom 
waveguides, for which the idealized limits a — > and 
N —>■ oo do not strictly apply a crossover must oc- 
cur from a effectively ID system, applicable when the 
waveguide is so long and narrow (high transverse fre- 
quency) that transverse excitations are frozen, to a 3D 



system with BEC accurately treated by the GP equa- 
tion (weaker transverse binding), the basis of most the- 
oretical work on trapped BECs. Detailed analysis of 
this crossover is important for comparison with exper- 
iments, since the ID regime has already been achieved 
experimentally |^|, |?], [|, |i| and the Tonks regime (ID 
and sufficiently large <?id) is being approached The 
dynamical reduction from 3D to ID and precise condi- 
tions on parameters necessary for achievement of both 
the ID limit and the Tonks-gas limit of the ID regime 
have been discussed in detail by Olshanii [10] and by 
Petrov et al. p| . In addition, Dunjko et al. [12| have in- 
vestigated the crossover between the Thomas-Fermi and 
Tonks-Girardeau regimes in a ID trap. Here we note 
only that the ID regime occurs when the waveguide is 
so thin (transverse frequency so high) and density and 
temperature so low that the longitudinal thermal and 
zero-point energies are both low compared with the low- 
est transverse excitation energy, resulting in "freezing 
out" all transverse excitations. Achievement of the Tonks 
limit requires, in addition, that the scattering length a 
is large enough and/or ID density n low enough that 
h 2 n/mgiD <C 1 where the effective ID coupling constant 
is giD — 2h 2 a/m£Q and £q = \fh/muJo is the transverse 
oscillator length. 

Trap geometry: A particularly convenient geometry for 
discussing the crossover is a toroidal trap of high aspect 
ratio R — L/£ . The transverse trap potential is sym- 
metric about an axis consisting of a circle on which the 
trap potential is minimum, and harmonic with respect to 
a coordinate p measured transversely with respect to this 
circle. This geometry can equally well be interpreted as 
an infinitely long, straight cylindrical waveguide with pe- 
riodic boundary conditions in the longitudinal direction. 
Toroidal traps of this form have been experimentally pro- 
duced and loaded |l| @. 

Hamiltonian: We use a many-body Hamiltonian with 
harmonic transverse binding and the usual Fermi pseu- 
dopotential interaction v(r i:) ■) = 47ra<5(r,y) with a positive 
s-wave scattering length a. This leads to a well-defined 
problem in ID, the LL model Q. Our toroidal system 
is "almost ID" since the transverse dimensions are con- 
fined, and we find that the variational problem with the 
Fermi pseudopotential does not encounter the difficul- 
ties (divergences and poorly-posed variational problem) 
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L5, 16| found in the 3D case. The Hamiltonian is then 
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where g$D = ^irafi 2 /m is the 3D coupling constant. The 
Laplacian is to be expressed in cylindrical coordinates 
r j = ( z j'Pji®j) wner e Zj, with < zj < L, is a ID co- 
ordinate measured around the torus circumference, pj is 
a transverse radial coordinate measured from the central 
circular torus axis, and 9j is the azimuthal angle about 
this axis. 

Variational ground state: We use a trial variational 
ground state which assumes factorization of longitudinal 
and transverse parts, with the transverse part depend- 
ing on a single transverse orbital (f> tr independent of az- 
imuthal angle: 

N 

$ (ri,---i r Jv) = ®long{zir-',ZN) Y[<Ptr{Pj) ■ (2) 

Use of a single transverse orbital is justified in two differ- 
ent limits: (a) tight transverse confinement, transverse 
excitations frozen, <p tr is the unperturbed transverse os- 
cillator ground state; (b) weak transverse confinement 
and low density, 4>tr is the transverse part of GP or- 
bital, assuming factorization of longitudinal and trans- 
verse parts of this GP orbital. Note, in connection with 
case (b), that a GP orbital well approximated by a Gaus- 
sian with respect to both z and p factorizes automati- 
cally, as does any Gaussian. More generally, we leave the 
functional form of cf>t r free, to be determined as part of 
the minimization of the variational ground state energy 
Eq = ($>o\H\&o). Assuming </> tr and $/ orig normalized 
according to 
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one finds that the energy expectation value of $o decom- 
poses as Eq — Netr + Ei ong with 

e tr = ^ 2npdp <f> tr [^- p g~ p Pg- p + ~^qP ) <t>tr 

Elo„g = / dZ\ ■ J dzN^ on gHi on g^i on g (4) 
JO JQ 

where Hi ong is an effective longitudinal Hamiltonian 
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and g\D is an effective ID coupling constant 

/>oo 

giD[M=93D 2lTpdp \<t>tr(p)\ 4 ■ (6) 

Jo 



Note that the separation of transverse and longitudinal 
energies is only partial, since gm depends on <j) tr . We 
take the value of the gm corresponding to the transverse 
ground state, with oscillator length t$, as a reference 
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and define the fractional ID coupling constant g\o — 
gm/g- For a ID system gn> = 1 whereas it decreases as 
the system crosses over to 3D. The relevant dimensionless 
intensive variable of a ID system 7 = (rngm) / '(h 1 'n) j^j 
suggests a dimensionless measure n = (7% 2 n)/(mg) of the 
linear density n = N/L. We should keep in mind that g oc 
<x>o and will vary as the transverse confinement changes. 

The total energy Eq is to be minimized with respect to 
variation of both <ptr and $; D ng subject to the normaliza- 
tion constraints. We imagine this done in two steps, first 
holding <ptr constant and minimizing Ei ong with respect 
to $/ ong (zi, • • • , zjv), then minimizing the resultant Eq 
with respect to <f> tr (p). For fixed 4> tr , hence fixed g\D, the 
global minimum of Ei ong is realized by the exact ground 
state of Hiong , which is the well-known LL Bethe Ansatz 
solution [Q. One may instead use some simpler varia- 
tional trial state for <&i ong , obtaining an upper bound to 
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to be minimized with respect to variation of <ptr subject 
to the normalization constraint: 
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where the transverse chemical potential ptr is the La- 
grange multiplier for the transverse normalization con- 
straint. Evaluation of the functional derivatives leads to 
the following generalized transverse GP equation: 
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where \itr is to be adjusted so that <f>t r satisfies the 
normalization constraint. The solution depends on g\o 
which in turn depends on 4> tr , so the solution has to be 
determined by a self-consistent iterative procedure, by 
making an initial guess for g\D, evaluating dei ong / dgm 
at this value of giD, solving Eq. (8) for c/) tr , determining 
a new gio from Eq. (6), and iterating to convergence. 

Two of us have previously developed a variational the- 
ory [ljj ] assuming the same toroidal geometry and using a 
variational trial state for $/ OTlg based on the variational 
pair theory of many-boson systems In that case 

the transverse GP equation (8) reduces to the previous 
one, Eq. (16) of ]Ti| , after correction of a typo therein 
p9| . That work was devoted to investigation of the BEC- 
Tonks crossover in the ID regime where the transverse 



orbital is frozen in the unperturbed transverse oscillator 
state, and no transverse GP equation was solved. Here 
we are concerned with a different crossover, namely the 
fD-3D crossover, which depends crucially on the solu- 
tion of the transverse GP equation. In the following two 
sections we shall separately consider first the case where 
*&ionci is approximated by the GA variational pair the- 
ory Q as in jl7| , which gives accurate results except in 
the Tonks gas regime (ID limit and very large scatter- 
ing length). Then we shall work out the solution using 
the LL theory |2| , the exact ground state of Hi ong ■ This 
latter is accurate even in the Tonks regime, but is more 
complicated to work out since the LL energy per par- 
ticle cll is known only from numerical solutions of the 
nonlinear LL integral equation. 

Pair theory solution: Approximating the longitudinal 
energy by the ID GA pair theory energy, we will first 
set eionn = ep, the expression for which we derived in 
Rcf. mil and can be written as 
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Here KAl — A) and E(l — A) are complete elliptic inte- 
grals Q an d / is the Bose-condensed fraction related to 
the dimensionless pair theory parameter A through the 
coupled equations 
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The partial derivative of the pair theory energy has to 
be taken before minimization [ ft7| rather than that of the 
expression in Eq. ( |Io| ) which is true only at the minimum 

J^ = f [2-/ 2 -2A/ 2 + A 2 / 2 ]. (12) 

LL solution: Next we shall use the LL energy cll 
which is the global minimum of Ei ong with respect to unre- 
stricted functional variation of <&i on g- It yields, together 
with Eq. (8), the best possible variational solution ob- 
tainable with a trial state of form (2): 
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where e("f) is obtained by solving the Lieb-Liniger system 
of integral equations H . The partial derivative 



LL 



dgiD 



2m hi 



e'( 7 ) 



(14) 



is the same regardless of the order in which the expecta- 
tion and the derivative are evaluated, unlike in pair the- 
ory since the LL solution is an exact eigenstate of Hi ong 
and thus obeys the Hellmann-Feynman theorem. 
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(10) FIG. 1: Left axis: ratios ell/ef and gm/g- Horizontal axis: 
logarithm of scaled linear density n = (h 2 n)/(mg). Right 



axis: logarithms of scaled energies calculated according to pair 
theory h 2 ep/(mg 2 ) and Lieb-Liniger theory h 2 ELL/(mg 2 ). 
The plots are shown for three different strengths of the trans- 
verse trap frequencies u>o directly proportional to the value of 
the dimensionless quantity mg 2 /(h 3 ujo) = (a) 10 -2 , (b) 10 2 
and (c) 10 8 . Note that the scale of actual linear densities n 
differ in the three plots since g is different for each plot. 

Numerical results and discussion: For both pair theory 
and Lieb-Liniger theory, the transverse GP equation was 
solved numerically using the method of discrete variable 
representation |2l| , a method based on Gauss quadrature 
with a well-defined discrete normalization. Iterating the 
solution of the GP equation alternately with the evalu- 
ation of the longitudinal energies lead to self-consistent 
values for guj and ti ong in both theories. In using LL 
theory we used the tabulated values of e(7) in |2^] for 
intermediate values of 7 and the limiting expressions in 
Q for high and low values of 7. For low 7, i.e. for high 
linear densities, the energies in both pair theory and LL 
theory coincide with the Bogoliubov energy ngm/2. 

In the Tonks limit of completely impenetrable bosons, 
all of the atoms are "fermionized" and the energy per 
atom is simply the energy of free fermions 
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We can get a sense of the degree of "fermionization" of 
the system by evaluating the ratio ell/^f- In the Tonks 
limit this would be unity, but as the atoms become pene- 
trable the energy approaches a linear dependence on the 
density so that the ratio will approach zero. In Fig. [I] we 
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transverse trap strength, Logjmg /fi co ) 

FIG. 2: Log-log plots of scaled linear densities h 2 n/mg ver- 
sus transverse confinement strength mg 2 /(h 3 u}o) at the BEC- 
Tonks crossover points defined by ei ong /eF = 0.5 and at the 
fD-3D crossover points defined by gio/g = 0.5 



plot this ratio (for LL theory) as a function of the den- 
sity n for different values of the dimensionless quantity 
mg 2 1 (h 3 luq) proportional to the transverse trapping fre- 
quency uj{). In the same figure we also plot gm = gm/g 
as function of n. While the plot of ei on g/eF gives a mea- 
sure of the impenetrability of the atoms, the plot of gm 
provides a measure of the dimensionality, since gio = 1 in 
ID and decreases as the system crosses over to 3D. Since 
g itself depends linearly on ujq the actual magnitudes of 
the densities are different in the three plots shown in 
Fig. 0. Along the right axis in Fig. |l| we plot the scaled 
longitudinal energies computed from pair theory and LL 
theory. We see that the density range where the pair 
theory energies breaks away from the LL energies cor- 
responds closely to the region where eiong/^F starts to 
deviate from unity. This is as we would expect since pair 
theory cannot describe the Tonks regime. 

We see that for low values of the transverse trapping 
potential as in Fig. |l](a) there is a regime of density where 
the gas is one dimensional but not yet a Tonks gas. Pair 
theory is valid in this region and it intrinsically allows 
for a "Bose-condensed" fraction / which has the more 
general interpretation as the non-trivial fraction of atoms 
in the ground state. Thus in this regime we should see a 
ID condensate, in the sense defined earlier. 

However, as we increase the transverse trapping poten- 
tial we see that the region where such a ID condensate 
can exist gets narrower as the curve for gio approaches 
the curve for ei on g/ep until the former is on the left of 
the latter as we see in Fig. 0(c); this means that for 
large transverse trapping potentials there will be never 
be a ID condensate and the system will pass directly 
from 3D to the impenetrable Tonks gas regime. This is 
clearly illustrated in Fig. ||; using guj = 0.5 as the cri- 
terion for crossover from ID to 3D and ei on g/^F = 0.5 
as the criterion for crossover to the Tonks gas, we plot 
the scaled densities at the crossover points as a func- 
tion of the transverse confinement strength measured 



by mg 2 / {h? uj ) . The triangular region between the two 
curves on the left of their point of intersection roughly 
defines the range of scaled densities and transverse trap 
frequencies over which a ID condensate can exist. On 
the right of the point of intersection the trap strength is 
too high. 

In conclusion, we have used Lieb-Liniger theory as well 
as pair theory to study the crossover of an axially homo- 
geneous 3D Bose condensed system to effective ID and 
eventually to an impenetrable regime. We have evalu- 
ated the densities at which the system crosses over from 
3D to ID and then to a Tonks gas for different trans- 
verse trap frequencies. We have demonstrated that for 
weak transverse confinement there is a physically allowed 
regime where a ID axially homogeneous condensate can 
exist while for sufficiently high transverse confinement 
the only possiblities are a 3D condensate or a Tonks gas. 

KKD thanks T. Bergeman for useful discussions. This 
work was supported by Office of Naval Research grant 
N00014-99-1-0806 and by the US Army Research office. 
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